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ABSTRACT 

Motivated by recent suggestions that a number of observed galaxy clusters have masses 
which are too high for their given redshift to occur naturally in a standard model 
cosmology, we use Extreme Value Statistics to construct confidence regions in the 
mass-redshift plane for the most extreme objects expected in the universe. We show 
how such a diagram not only provides a way of potentially ruling out the concordance 
cosmology, but also allows us to differentiate between alternative models of enhanced 
structure formation. We compare our theoretical prediction with observations, placing 
currently observed high and low redshift clusters on a mass-redshift diagram and find 
- provided we consider the full sky to avoid a posteriori selection effects - that none 
are in significant tension with concordance cosmology. 

Key words: methods: analytical - methods: statistical - dark matter - large-scale 
structure of Universe - galaxies: clusters 
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1 INTRODUCTION 

In the standard "concordance" (ACDM) model of cos- 
mology, structure formation proceeds in a hierarchical, 
'bottom-up' fashion, with small scale perturbations in the 
initial distribution of Cold Dark Matter (CDM) collapsing 
first, before merging over time to form larger and larger 
haloes. The initial, nearly Gaussian, distribution of these 
overdensities means that large fluctuations, and correspond- 
ingly large halo masses, should be extremely rare in the 
early universe. However, many of the plausible extensions 
to the concordance model have been found to be capable 
of enhancing (or depleting) st ructure formation, includ- 
ing primordial non-G aussianity dLucchin fc Matarresel[l988| ; 
Pillcpich ct al. 2010), modified gravity JSchmidt et ah 2009; 
Ferraro et alj l201ll) and scalar field dBaldi fc Pettorinol 
201ll ; lMortonson et alj|201ll ; [Tarrant et alj|201ll ) scenarios. 



Furthermore, each of these extensions can affect structure 
formation in different ways at different times in the history 
of the universe, meaning that if we can construct a history 
of the growth of structure we can discriminate between com- 
peting models. 

In addition, much attention has recently been paid 
to the prospect that the existence of a single extreme 
(in terms of its high mass and redshift) object in the 
universe has the potential to rule out at a high sig- 
nificance level cosmological models in which it is cor- 
respondingly unlikely to exist. Following observations of 
a series of apparen tl y extreme object s dJee et al. 20091; 



a series 01 apparen tl y extreme object s (IJee et al.l 2009; 
iBrodwin etail l2010t iFolev et all 1 20 lit iMenanteau et all 



l201ll ; iPlanck Collaboration et al.ll201ll ; ISantos et al]|201ll ). 
some authors have claimed t hat such object s are 
either highly unlikely t o exist (Ijimenez fc Verdd 120091; 
Holz fc Perlmutterl |2010| ; ICavon et al.l |201ll ; iHovle et al.l 



20 111 ) in a concordance ACDM cosmology or, more power- 
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fully, are significantly larger th an the expected most massiv e 
object in a ACDM unive r se JChongchitnan fc Silkl l201ll l. 
although Waiz mann et al.l (|2011aT ) contend this result. As 
well as highlighting tensions with the standard model, 
various authors have shown how such high-mass, high- 
redshift objects may be explained by an enhanced rate of 
structure formation generated by the inclusion of primor- 
dial non-Gaussianity (jEnqvist et alJl201ltlHovle et al.ll201ll ; 
ICavon et alJl201ltlChongchitnan fc Silkll201ll) or scalar field 
dark energy ( Waizmann et all 2011bl ) into the concordance 
cosmology. 

In this paper we b uild on our previous wo rk using ex- 
treme value statistics (|Harrison fc Coles! 120111 ) to test the 
current concordance model with extreme galaxy clusters 
as well investigating the different behaviours of extreme 
clusters in alternative models. We construct the probabil- 
ity contours of the highest mass cluster for all redshifts in 
an observational survey and directly compare the current 
best fit ACDM model with observations showing that, if we 
consider the full sky (in order to avoid a posteriori selection 
effects) then none of the currently observed extreme galaxy 
clusters are in tension with the ACDM mod el. Then, using 
information from the CoDECS (|Bardill2011ah N-body simu- 
lations, we show how, if any future observations do exclude 
the ACDM model, extreme clusters can potentially be used 
to understand which of the alternative models available may 
best explain the enhanced structure formation. 



© 2012 RAS 



2 Harrison & Coles 



The paper is organised as follows. In section [2] we in- 
troduce extreme value statistics and show how they may be 
used to predict the mass of the most extreme cluster within 
an observational survey. Section [3] compares predictions for 
a ACDM universe with observations. Section [4] shows how 
the plot of extreme value contours against redshift can high- 
light differences between cosmological models. In section [5] 
we conclude and discuss prospects for future work in this 
area. 



2 METHOD 



Extreme Value Stati stics (EVS) (|Gumbell 1 19581 ; 
iKatz fc Nadaraiabl 120021 ') seeks to make predictions for 
the properties of the greatest (or least) valued random 
variable drawn from an underlying distribution. If we 
consider a sequence of N random variates {Mi} drawn 
from a cumulative distribution F(m) then there will be a 
largest value of the sequence: M max = sup{Mi, . . . Mjv}. If 
these variables are mutually independent and identically 
distributed then the probability that all of the deviates are 
less than or equal to some m is given by: 



$(M ma * ^m;N) 



= F 1 (M 1 ^m)...F N {M N 
= F N (m) 



< m) 



and the probability density function (pdf) for M n 
found by differentiating fT]): 



4>(M B 



m;N) = Nf(m) [F(m)Y 



(1) 
is then 

(2) 



where f(m) = dF(m)/dm is the pdf of the original distri- 
bution. This gives the exact extreme value pdf for N ob- 
servations drawn from a known distribution F(m); for more 
discussion of the advantages of this approach over others in- 
volving the use of asymptotic theory, see 
citeHarrison2011. To apply this general result to the con- 
crete example of the most massive cluster we use the appro- 
priate halo mass function, which gives the number density 
of haloes n(M), to derive f(m) and F(m) according to: 



f(m) 
F(m) 



dm 



1 dn(m) 
ntot 
1 

ntot 



dM 



dn(M) 



dM 



where the normalisation factor 
dn{M) 



ntot 



dM- 



dM 



(3) 
(4) 

(5) 



is the total (co-moving) number density of haloes. For a 
constant redshift box of volume V the total number of ex- 
pected haloes N is given by n to tV . These distributions can 
be inserted into equation ([2)| to predict the pdf of the highest 
mass dark matter halo within the volume. Although this pro- 
cedure expli citly assumes that halo es are uncorrelated, we 
have found ([Harrison fc Colesj|201ll ) that the results closely 
match those of lDavis et al.l ( 20111 ). who construct their ex- 
treme value distribution as the differential of the void prob- 
ability: 



dP (m) 
dm 



(6) 



(allowing them to account for complications due to halo cor- 
relations and bias), for the high masses of relevance to the 
inference we are interested in. 

In a cosmological survey we observe clusters at vari- 
ous redshifts along our past light cone rather than on a 
single spatial hypersurface at fixed z. If we wish to con- 
struct the EVS for galaxy clusters within an observational 
survey which covers a fraction / a ky of the sky between red- 
shifts 2 m i n and 2 max we therefore need to take into account 
both the effect of the growth of structure with decreasing 
redshift on the halo mass function n(m, z) and the observa- 
tional volume we are probing in an expanding universe, via 
the volume element dV/dz. Doing this allows us to form the 
pdf of halo masses within that survey as: 



f(m) 


/sky 

Ntot 


F(m) 


/sky 

Ntot 


where 






= /sky J 



dz 



dV dn(m, z) 



dz dm 



dz dM 



dV dn(M, z) 



dz dM 



dz dM 



dV dn(M, z) 



dz dM 



(7) 
(8) 

(9) 



and then feed these distributions into our extreme value 
prescription (T2J) (of course it is impractical to integrate nu- 
merically to infinite endpoints and so finite limits of 12 < 
log 10 m < 18 are chosen; we have checked that this choice 
makes no difference to the conclusions). In order to make 
best use of this information, we want to be able to see the 
distributions for all redshifts at once; we hence construct the 
EVS distribution for narrow bins in redshift space Az = 0.02 
(chosen so that Nbms » N c i U3t ers and the highest expected 
mass for all redshifts remains the same as for Nbins = 1), 
integrate over these pdfs to find the 66%, 95% and 99% con- 
fidence regions and plot these, along with the peak of the 
distribution, for all redshifts < z < 2. This can then be 
used to test the cosmological model: if an observed cluster 
lies above e.g. the 95% region of such a distribution, then 
we may say we have a correspondingly significant detection 
of enhanced structure formation at that redshift. 



3 COMPARISON WITH OBSERVATIONS 

We can now apply this technique to find out if any cur- 
rently observed clusters are discordant with the concordance 
ACDM predictions. We emphasize that, because we are pre- 
dicting the distributions of the most massive cluster at each 
redshift, if even a single galaxy cluster lying outside the ex- 
treme value contours when placed on a mass-redshift plot 
can be seen as a significant detection of deviation from con- 
cordance cosmology. 



3.1 Calibration of Cosmology and Cluster Masses 

In order to meaningfully compare our theoretical predic- 
tions to observations, we need to carefully ensure our ACDM 
cosmology is correctly calibrated. As the ingredients for 
our concordance cosmology here, we use a linear mat- 
ter power spectrum P(k) calculated using the numerical 
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Cluster 



A/fEdd 
iu 200m 



/M 



x 10 



A2163 1 


0.203 


3.04 


+0.87 

-0.67 


X 


10 15 


A370 1 


0.375 


2.62 


+0.87 
-0.67 


X 


10 15 


RXJ1347 1 


0.451 


2.14 


+0.60 
-0.48 


X 


10 15 


ACT-CL J0102 2 


0.87 


1.85 


+0.42 
-0.33 


X 


10 15 


PLCK G266 3 


0.94 


1.45 


+0.27 
-0.20 


X 


10 15 


SPT-CL J2106 4 


1.132 


1.11 


+0.24 
-0.20 


X 


10 15 


SPT-CL J0546 5 


1.067 


7.80 


+ 1.27 
-0.90 


X 


10 14 


XXMU J2235 6 


1.4 


6.82 


+ 1.52 

-1.23 


X 


10 14 


XXMU J0044 7 


1.579 


4.02 


+ 0.88 
-0.73 


X 


10 14 



Table 1 ■ The extreme clus ters c onsidered in th i s pa- 
p er ([11 iMaughan et alj l201l | , [ 21 iMenanteau et ail l201ll. 
[3lPlanck Collaboration et al.l l201ll, [41 iFolev et ail l201ll, [5 1 
iBrodwin et ail l20ld , [61 I Jee et al.1 12009| . [71 ISantos et alj|201ll). 

A^JaOm is calculated using the numerical code of l|Zhao et ail 
2009) to convert from M20OC (where necessary) and equation i ] 1 12 D 
to include the Eddington bias. 



code CA Mffl and the WMAP 7+BAO+H0 Mean paramet- 
ers from iKomatsu et ail |201ll ). From this we calculate the 
variance of the matter field, smoothed with a top hat win- 
dow W(k; R) of radius R = (3m/47rp„ li o) 1 '' 3 , evolved to a 
redshift z with the linear growth function D+(z) (normal- 
ised to D+(0) = 1): 

/•oo flu 

o 2 (m,z) = D 2 + (z) jp k 2 P(k)W 2 (k; R). (10) 
Jo 2tt 

This is used as the input for the halo mass function from 
iTinker et ail |2008l ): 



dn(m, 
dm 



= A 



+ 1 



c/a 2 pm,o dln(a 1 ) 
m dm 



(11) 



where p mi o is the mean density in the universe at red- 
shift z = 0, 8c — 1.686 is the critical overdensity for col- 
lapse and {A, a, b, c} are given the values for m,200m (the 
mass corresponding to the portion o f the cluster which has 
density greater than 200p m ,o) from ITinker et ail (|2008l ) of 
{0.186,1.47,2.57,1.19}. 

When comparing to real-world clusters we need to cor- 
rect for the fact that theoretical mass functions are defined 
with respect to the average matter density p m ,z, but observ- 
ers frequently report cluster masses with respect to the crit- 
ic al density p c - In o r der to do th is, we follow the procedu res 
of lWaizmann et"ail (|2011al ) and iMortonson et ail (|201ll ) to 
convert all cluster masses to m,2oom , and correct for Edding- 
ton bias. Eddington bias refers to the fact that there is a 
larger population of small mass haloes which may up-scatter 
into our observations than there are high mass haloes which 
may down-scatter into them, and is corrected for using: 



In m 



1 _l 1 2 

2 e<7lnm ' 



(12) 



where e is the local slope of the halo mass function and of n m 
is the measurement uncertainty for the cluster mass. 

In order to ensure we are avoiding a posteriori selection 
(by only performing our test in regions in which we have 
already observed something which we believe to be unusual) 
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Figure 1. Extreme value contours and modal highest-mass 
cluster with redshift for a ACDM cosmology, along with a set 
of currently observed 'extreme' galaxy clusters. None lie in the 
region above the 99% contour and hence are consistent with 

a concordance cosmology. 

we set f a ky = 1. This is both the most conservative estimate 
and, we believe, the correct one for testing 'the most extreme 
clusters in the sky'. 



3.2 Results 

We now seek to use the apparatus described above to 
test if any currently observed objects are significantly ex- 
treme to give us cause to question ACDM cosmology. We 
consider the set of recently observed, potentially extreme 
clusters shown in table 1 in a ACDM cosmology as described 
above. The extreme value contours (light - 99%, medium - 
95%, dark - 66%), most likely maximum mass M^ ax (solid 
line) and the cluster masses and redshifts (stars) are plot- 
ted in Figure [T] The plot shows the expected features of 
a peak in maximum halo mass at z ~ 0.2 (the location 
and hei ght of which is in broad agreement with the ana- 
lysis of IHoIz fc Perlmutterll2(5l0l ). As can be seen, none of 
the currently observed clusters lie outside the 99% confid- 
ence regions of the plot meaning that there is no current 
strong evidence for a need to modify the ACDM concordance 
model from high-mass high-redshift clust ers. This appears 
to be in agreement with the findings of IWaizmann et al] 
il2011aD for a similar se t of cl usters, but in contradiction 
to IChongchitnan fc Silkl l|201ll fl who find that the cluster 
XMMUJ0044 is a 4a result (i.e. should lie well above the 
99% region in Figure [TJ , whilst here we find it to be well 
within the acceptable region. 



4 TESTING COSMOLOGICAL MODELS WITH 
EXTREME CLUSTERS 

In addition to simply ruling out ACDM cosmology with 
massive clusters, we may also consider whether extreme ob- 
jects offer the potential to discriminate between different 



2 In an updated version of this analysis, Chongchitnan 
find no tension with ACDM 



Silk 
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alternative models. Whilst many alternative models are cap- 
able of predicting enhanced structure formation, the exact 
scale and time dependence of the enhancement will differ 
from model to model. Here we consider two models which 
have a well defined and investigated effect on the halo mass 
function, and hence are relatively simple to calculate the 
extreme value statistics over a range of redshift for: local 
form primordial non-Gaussianity and the bouncing, coupled 
scalar field d ark energy model labelled as 'SUGRA003' in 
iBaldil |2011bh . These should be regarded as toy models - 
our aim is to show how the extreme value statistics can be 
used to select between different models, rather than make 
definite predictions. 



4.1 Models Considered 

We make use of the CoDECS s i mulati on s kind l y mad e 
publicly available by iBaldi et al.l |2010r i; IBaldil ^Ollal ). 
This suite of large N-body simulations includes realisa- 
tions of both ACDM and a number of coupled dark en- 
ergy cosmologies. Here, we compare the CoDECS ACDM- 
L (where 'L' is for 'Large') simulation of the concord- 
ance cosmology to both the primordial non-Gaussianity 
and the SUGRA003 (supergravity) bouncing dark energy 
models. Primordial non-Gaussianity, motivated by consid- 
erations of the fluctuations of the inflaton field, is one of 
the most widely exp lored modifications to the concord- 
ance cosmology (e.g. iDesiacques fc Seliakl I2010T I and has 
long (|Lucchin fe Matarresd Il988f ) been known to affect the 
abundances of high-mass gal axy clusters. It has also been the 
model most invoked (Ijirnenez fe Verdd 20091; ICavon et alj 



1201 ll : IChongchitnan fe Silkll201ll ; iHovle et ail 1201 lh to~acT 
count for apparently over-massive high redshift objects, all 
of these authors reporting values of local non-Gaussianity 
parameter /jvl ~ 300 — 500 as being able to account for 

such clusters. 

However, IBaldil (|2011bh points out that such models 
enhance numbers of high mass clusters at all redshifts, cre- 
ating tension with observations at low redshift in the at- 
tempt to alleviate them at high redshift. As an alternative 
scenario, the supergra vity-motivated scalar field scenario of 
iBrax fe Martini (| 19991 ) is considered. This model includes a 
scalar field component <j> which couples to dark matter with 
a coupling strength j3 and has the self interacting potential: 



v{<(>) 



-a <A z /2 



(13) 



This scalar field component acts as a 'bouncing' dark en- 
ergy; structure formation is enhanced at early times, but 
is suppressed with respect to ACDM after the point at 
which the evolution of (f> changes sign (the 'bounce'), mean- 
ing ACDM values for os can still be reproduced at z = 
0. In order to mat ch background obser vables given by 
WMAP7 constraints (|Komatsu et al.ll201lf ). the SUGRA003 
version of the potential has 'tuned' parameters {B,a,/3} = 
{0.0202,2.15,-0.15}. 

For ACDM and SUGRA mode ls, we fit a halo mass 
function of the I Tinker et al.l |2008l 'l form directly to the 
haloes identified using a friends-of-friends (FoF) algorithm 
with linking length I = 0.2d (where d is the mean inter- 
particle separation) in the relevant CoDECS simulation 
(ACDM-L and SUGRA003-L respectively). For the primor- 



dial non-Gaussianity model, we apply a non-Gaussian cor- 
rection factor 7?.(/nl) to the halo mass functi ons found in the 
ACDM -L simulation, choosing the 7?.(/nl) of lLo Verde et ail 
(|2008h (LMSV): 



-1 + 



dS 3 



- 1 



(14) 



TZ-LMSvifNL) = 

where S3 is the normalised skewness of the matter density 
field, for which we use the approximation: 

S 3 ~ 3 x lO^fNLCj- 1 (15) 

given by equation (2.7) of lEnovist et al.l (|201ll ). We adopt a 
value of /jvl = 300 for our non-Gaussian model as it is both 
consistent with the observational findings discussed above 
and leads to a similar magnitude of enhancement of struc- 
ture formation at high redshifts as the SUGRA003 model. 

The values of H{z) and D+(z) required to find dV/dz 
for all three models are calculated using the tabulated 
growth functions and expansion histories for the cosmolo- 
gies, numerically calculated from the evolution equations 
and provided on the CoDECS website. 

4.2 Results 

With the halo mass functions and expansion histories of 
each cosmology we are then able to carry out the procedure 
of section [5] to find the EVS of objects within an obser- 
vational survey in each cosmology, the results of which are 
shown in Figure [2] Plotted are extreme value contours (light 
- 99%, medium - 95%, dark - 66%) for the ACDM model and 
the edges of the three extreme value contours for the non- 
Gaussian and SUGRA models (dashed lines) as well as the 
enhancement in the most likely maximum mass M^ ax over 
the ACDM predictions. As can be seen (and as expected) 
the primordial non-Gaussianity model shows an enhance- 
ment of the mass of the highest mass cluster at all redshifts, 
whilst the SUGRA model is capable of enhancing M^ ax at 
high redshifts whilst leaving it unchanged at more recent 
times. Thus, if ACDM is ruled out by both high and low 
redshift clusters, primordial non-Gaussianity could be seen 
as the favoured explanation whilst, if only high redshift ob- 
servations appear in contradiction, both non-Gaussian and 
SUGRA models would be allowed (unless the limit of an 
ideal, complete survey was reached). 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have presented a theoretical framework for 
the interpretation of extremely massive clusters in cosmo- 
logical surveys. We avoid ambiguities in the volume being 
probed by observed clusters which occur in the alternative 
approach of using single clusters to estimate the halo mass 
function n(m). Our method, in contrast, does not require 
such normalisation and is both more robust and, inevitably, 
more conservative, a combination of attributes which we be- 
lieve should best describe confrontations of theory with ob- 
servations. 

We have used this method to show that existing ob- 
servations of extreme clusters are not incompatible with 
the predictions of ACDM concordance cosmology. However 
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Figure 2. Extreme value contours for ACDM (shaded regions), 
/jvl and SUGRA models (dashed lines). Lower plots are the en- 
hancement of modal highest-mass cluster over the ACDM value, 
showing different behaviour for the two alternative models. 

in anticipation of future, more complete surveys, we have 
also presented a definitive calculation of the region in the 
mass-redshift plane in which the presence of a single cluster 
would exclude the standard model with high confidence. We 
have also shown that different cosmological models includ- 
ing physics beyond the concordance picture are capable of 
making significantly different predictions of the allowable 
regions in the mass-redshift plane. This means that whilst 
current extreme cluster data do not discriminate between 
the ACDM model and these alternatives, future data may 
well be capable of doing so. We intend in due course to ex- 
tend the approach presented here by incorporating it into a 
full model selection analysis. 
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